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States of Superconductivity and Magnetism 
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Ginzburg-Landau (GL) equations for the coexistent states of superconductivity and mag- 
netism are derived microscopically from the extended Hubbard model with on-site repulsive 
and nearest-neighbor attractive interactions. In the derived GL free energy a cubic term that 
couples the spin-singlet and spin-triplet components of superconducting order parameters 
(SCOP) with magnetization exists. This term gives rise to a spin-triplet SCOP near the inter- 
face between a spin-singlet superconductor and a ferromagnet, consistent with previous the- 
oretical studies based on the Bogoliubov de Gennes method and the quasiclassical Green's 
function theory. In coexistent states of singlet superconductivity and antiferromagnetism it 
leads to the occurrence of n- triplet SCOPs. 
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1. Introduction 

The coexistence and competition of superconductivity and magnetism have been impor- 
tant issues in various strongly correlated electron systems, e.g., high-T c cuprate supercon- 
ductors. 1} This is because these two ordered states originate from the same interaction; thus, 
understanding their relation may give insight into the mechanism of superconductivity. 

Heterostructures composed of superconductors and magnetic materials may be useful 
systems for studying these problems. The properties of the states near an interface strongly 
depend on the materials used, especially the symmetry of superconducting (SC) states and 
the underlying electronic states, 2 ~ 16) namely, the shape of Fermi surfaces and the type of in- 
teractions. The interface states of heterostructures may have quite different characters from 
those in the bulk. Not only the coexistence of the original order parameters (OPs), but also 
new ordered states may arise depending on the constituent material. For example, spin-triplet 
SCOPs can occur near the interface between a spin-singlet superconductor and a ferromag- 
net. This was theoretically found using the Bogoliubov de Gennes (BdG) method 5) and the 
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quasiclassical Green's functions theory. 9 ' 10) Therefore, by combining various types of super- 
conductors with ferromagnets or antiferromagnets, we may know the conditions under which 
a particular SC state can be realized. 

In this paper, we derive GL equations and the GL free energy microscopically from a 
tight-binding model on a square lattice with on-site repulsive and nearest-neighbor attractive 
interactions, i.e., the extended Hubbard model. Although this model is a minimal one for 
treating magnetism and unconventional superconductivity, it exhibits s-, d-, and chiral (p x ± 
/p y )-wave superconductivity, 17 ' 18) and ferro- and anti-ferromagnetism 19) for different choices 
of the parameters, especially the electron density (in other words, the shape of the Fermi 
surface). For this reason this model may be used to examine the material dependence of 
the interface states of heterostructures composed of superconductors and magnetic materials. 
The method of deriving GL equations is based on that by Gor'kov 20) with an extension to 
include magnetic OPs. The resulting GL equations are coupled equations for all types of OPs 
including magnetization. (The GL equations for superconductors with s- and rf-wave SCOPs 
have already been obtained from a similar model. 21 ' 22) ) 

Although the GL theory is reliable only qualitatively except near T c , it can give a simple 
and clear description of the coexistence and competition of multiple OPs. Thus, it is comple- 
mentary to more sophisticated methods such as the BdG and quasiclassical Green's function 
theory. 

This paper is organized as follows. In §2, we present the model and treat it by mean- 
field approximation (MFA). In §3, GL equations and the GL free energy are derived for the 
coexistent states of superconductivity and ferromagnetism. The case of antiferromagnetism 
and superconductivity is examined in §4. Section 5 is devoted to summary and discussion. 

2. Model and Mean-Field Approximation 

We consider the extended Hubbard model on a square lattice, i.e., a tight-binding model 
that has on-site repulsive and nearest-neighbor attractive interactions. (We use the units fi = 
ks = 1, and the lattice constant is taken to be unity.) By treating the latter interaction using 
MFA, s-, d-, and chiral (p x + /p v )-wave SC states can be realized depending on the electron 
density. Namely, the symmetry of the SC state may be determined by the shape of the Fermi 
surface. 1718) Similarly, the magnetic order, either ferromagnetic and antiferromagnetic (AF), 
can be obtained by treating the repulsive interaction by MFA. 19) 
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The Hamiltonian of our model is given by 

H = ~t Yj E ^.U,r<>^ • + h. C ) - fio X C^C 



j j " 

where /i , U, and V are the transfer integral, chemical potential, and the on-site repulsive and 
nearest-neighbor attractive interactions, respectively; <x =T, J, is the spin index. The magnetic 
field is taken into account using the Peierls phase (pj C = ^ JT A • dl, with A and O = §f 
being the vector potential and flux quantum, respectively. We treat this Hamiltonian using the 
standard procedure of MFA: 

n n n }i -> (n n ) njl + {nj X )n n - (n n )(n n ), 
n n n a -» Aj,ec}/ n + A* f c yT c a - |A^| 2 . 
The SCOPs and magnetization (i.e., the OP for magnetism) are defined as 

A# = (cftc ei ), nij = -{n n - nji). (3) 

Then the mean-field Hamiltonian is written as 

Hmf E 'V^V^ " 



jcr 6=±x,±y 

j 

~ y Z E [A 7 j +5 c ; t +(5i c ; t T +/i.c.]+ J E , 

y <5= ±i,±> i 



(4) 



where 



i2 



(5) 



;' j <5=±.i,±y 

with n. } = (riff + n^) being the electron density at the site j. Here, ji = /j. - Un /2 is the 
renormalized chemical potential with n being the average electron density of the system. 

In order to derive the GL equations for OPs, we introduce the following thermal Green's 
functions: 

GAj, €, r) = -<r T c>(T)4>, F^,a €, r) = -(T^t)^,). (6) 

The equations of motion for G a and F^, (i.e., the Gor'kov equations) are obtained by taking 
their r derivatives and carrying out Fourier transformation to the Matsubara frequency is,, 
(= (2n + l)inT; T being the temperature). These equations can be transformed to the following 
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(7) 



coupled equations for and F^,: 

G T (j, e, is n ) = G Q (j, £, is n ) + G (j, k, ie n )A kMS F f u (k + 6, €, is„) 

k,S 

-U^j GoU, k, ie n )m fc G T (fc, t, is n ), 
k 

G±(j, £, ie n ) = G Q (J, €, is n ) G °(j' i£ n )A k+s ^ n {k + 6, 1, is n ) 

k,S 

+U G (j, k, is n )m k Gi(k, I, is n ), 

k 

F| T (j, €, is n ) = -V ^ G (k, j, -is n )A* k+6k G^(k + 6, £, is n ) 

k,S 

+u /] Go(k, j, -is n )m k F^{k, I, is n ), 
k 

F^ij, e, is n ) = V ^ G (k, j, -is^Al ^gG^k + 6, €, ie n ) 

k,S 

-U^j G (k, j, -is n )m k F f n (k, €, is n ), 

k 

where the summation on 6 (k) is over +x and ±y (all sites). Here, Go(j, t, ico n ) is Green's 
function for the system without A and m but with A satisfying 

(is„ + fx)G (j, t, ie n ) + t ^ Go(j + S, £, ie n )e^'^ = 6 U . (8) 

G is related to Green's function for the system without A, G , as G (j, €, is n ) = 
Go(j,€,iE n )e^j' e . Go(j,£,is n ) is the Fourier transform of G (p, is n ) = l/(is n - % p ) with 
% p = -2?(cos p x + cos Py) - pi. 

Spin-singlet and spin-triplet SCOPs on the bond (j, j+rj) are expressed in terms of Green's 
functions F^ and F^, 

(Af 0'))* = ^cnc Mi - c n c m r = \{A jtj+n + A j+tlJ y 

1 £ " 1 

(A ( P(j)T = ^{c n c Ml + c }i c m )* = -(A jJ+r] - Aj +r] j)* 

and the magnetization is similarly given using Gf and Gj as 

m j = |< C W ~ c )i c jO 

We substitute eq. (7) into eqs. (9) and (10) iteratively and keep the terms up to the third order 
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in OPs to get the following GL equations: 

(Af \j)Y = J] L(1) 0'' k > n, S)(Af\k))* + J] L {2 \j, k, k', r,, 6)(A ( P(k)Tm t 

k,S k.k',S 



/Uf 5,<5' ,5" 



-(Af (k))XAP(k'MA%)(k")y - (Af\k)nA^(k'))(Af\k'')T (1 1} 



+ 



+(Af ) (fe))*(Af(fe'))(A?(fe"))*] 
^ L (4) (j, *, fe", 77, S%Lf{k)ym w m w , 



k,k',k",6 



(A^\j)T = - E L(1) °'' 775 - E L(2) °'' *' 1> «5)( A ? ) (*))*»f 

+ E E L(3) °'' *• ' 775 5 ' ^")[(Af ) W)*(A^ ) (fc'))(A?(r ))* 



k,k',k" 6,6',6" 

-(Af \k)nA { "\k'))(A^(k")T - (A ( P(k)y(A { *\k'))(A%,\k'')y (12) 



+(Af ) (fc))*(A< 7 ; ) (fc'))(A^(fc"))*] 
^ L (4) (j, *, /V, 8)(Af\k))*m k >m k », 



k,k',k",S 



mj = Y J L (5 \j,k)m k + L^a^^.tf.^K^XAfCW-A^^XAgW] 
+ 2^ L a \j,k,k' ,k")m k m t m k „ 

Kk'.k" 

+ E E^^^^'^^'^K^^^^^)* - A^^xA^^or^, 



(13) 

where the functions L (n) (n = 1, • • • , 8) are given in Appendix A. 

From eqs. (11)-(13), it is seen that the equations for A (5) , A (r) , and m have the second- 
order terms of the forms mA (T \ mA (S \ and A (S) A (r) , respectively. This implies that the GL 
free energy should have the cubic term of the form mA (S) A (r) , and it is actually the case as we 
will see in the following sections. It should be noted that eqs. (11)-(13) are valid even when 
the OPs have rapid spatial variations, because we have not yet taken a continuum limit. This 
property is important when we consider the antiferromagnetic case in §4. 

3. GL Equations for Coexistent States of Superconductivity and Ferromagnetism 

In this section, we consider the coexistent states of superconductivity and ferromagnetism. 
The GL equation for the SCOP of each symmetry can be obtained by making a linear combi- 
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nation of eqs. (11) and (12): 

^U) = \ £ Af 0), A d (j) = \[ X AfCi) - E A^(j)], 

^=±i,±>i rj=±& rj=±y (14) 

A p ,C/) = i[A^C/)-A3(/)]. A w C/) = i[Aj n C/)-A5C/)]. 
Assuming that the SCOPs and magnetization are slowly varying, we take a continuum limit. 
The SCOPs in the linear and quadratic terms are expanded in powers of derivatives by denot- 
ing tj ^r,r k ^r': 

A d (k) -»A tf (r') 



A^r) + (r' - r)„V p A,(r) + -(r' - r) p (r' - r^V.A^r), 



1 . . (15) 

where the summations over jj. and v are assumed, and a similar approximation is carried out 
for m. The Peierls phase coming from G is also expanded in powers of A. Using the ap- 
proximation (p k j ~ — ^(r' — r) ■ A(r), the derivatives and A are combined to construct the 

2ni 

gauge-invariant gradient acting on A, D = V + — A, and we keep the terms up to the sec- 

00 

ond order in D. As a typical example for treating the derivative terms, the derivation of the 
second-order term for A s is presented in Appendix B. In the third-order terms, we neglect the 
derivative terms and the vector potential A as usual, namely, A (S) , A (r) , and m with the argu- 
ments k, k', and k" are replaced with A (5) (r), A (T) (r), and m(r), respectively. Rewriting A (S) 
(A (r) ) using A s and A d {A px and A py ), we carry out straightforward but lengthy calculations to 
get the following GL equations for SCOPs and magnetization: 



arA + 2A|A,.| 2 A, - K S (D 2 X + D 2 y )A s - K sd (D 2 x - D 2 y )A d 
+ K spm [(V x m) A px + (y y m)A py + 2m(D x A px + D y A py )] 
+ 7l | A d \ 2 A s + 2y 2 A 2 A: + y 3 (\A px \ 2 + \A py \ 2 )A s + 2y 5 (A 2 px + A 2 py )Al 
+ TvOApJ 2 - |A pv | 2 )A d + y,(A 2 px - A 2 py )A* d + y ms m 2 A s = 0, 



(16) 



a d A d + 2/3 d \A d \ 2 A d - K d (D 2 x + D 2 )A d - K sd {D 2 x - D 2 )A S 
+ K dpm [(V x m)A px - (V v m)A /n , + 2m(D x A px - D y A py )] 
+ 7l \A s \ 2 A d + 2y 2 A 2 A d + y 4 (\A px \ 2 + \A py \ 2 )A d + 2y 6 (A 2 px + A 2 py )A* a 
+ y 7 (\A px \ 2 - \A py \ 2 )A s + y s (A 2 px - A 2 py )A* s + y md m 2 A d = 0, 



(17) 
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a p A px(y) + 2{i p \A px (y)\ A px(y) - K pi D x(y) A px(y) - K p2 D y(x) A px(y) 

- (K p3 + K p4 )D x D y A py(x) - K spm [(V x(y) m)A s + 2mD x{y) A s ] 

- K dpm [(V x(y) m)A d + 2mD x(y) A d ] + y pl \A pyix) \ 2 A px(y) + 2y p2 A 2 py{x) A* px{y) (18) 
+ y 3 \A s \ 2 A px(y) + y 4 \A d \ 2 A px(y) + 2y 5 A 2 s A* px(y) + 2y 6 A 2 A* px(y) 

± y 7 (A s A* d + c.c.)A px(y) ± 2y^A s A d A* px(y) + y mp m 2 A px(y) = 0, 

a m m + 2/3 m m 3 - K m (V 2 x + V 2 )m 

+ l -K spm \A s {(D x A P xT + (DyApy)* } - {A* D X A S + A* D y A s ] + c.c.l 

\ (19) 
+ ^K dpm [A d {(D x A px y - (DyApy)*} - {A* px D x A d - A* py D y A d } + c.c] 

+ y ms m\A s \ 2 + y md m\A d \ 2 + y mp m(\A px \ 2 + \A py \ 2 ) = 0, 
where the coefficients appearing in eqs. (16)-(19) are given in Appendix C. 

Equations (16)-(19) are the coupled equations that determine the SCOPs and magneti- 
zation self-consistently. The most important point is that the second-order terms with a first- 
order derivative exist in the GL equations. They can induce triplet (singlet) SCOPs in a singlet 
(triplet) superconductor once the magnetization coexists inhomogeneously. It should also be 
noted that the coefficients in GL equations are determined microscopically, reflecting the na- 
ture of the electronic states of the original model, e.g., the shape of the Fermi surface. This 
property can be used to study the coexistent states of realistic materials to be considered. 

The GL free energy F up to the fourth order in OPs can be obtained from the above GL 
equations in such a way that the variations of F with respect to OPs reproduce eqs. (16)-(19). 
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The results are written as follows: 

F = F s + F T + F ST + F M + F SM + F TM + F STM , 

F s = J d 2 r[a s \A s \ 2 + /3 s \A s f + K s 0A s \ 2 + a d \A d \ 2 + /3 d \A d \ 4 + K d 0A d \ 2 

+ ri |A. s | 2 |A f/ | 2 +r 2 (A2(A:) 2 + c.c.) 
+K ds ((D x A<d(D x A s y - (p y A d )(D y A s y + c.c.)], 
F T = J d 2 r[a p (\ A px \ 2 + \A py \ 2 ) + f3 p (\A px \ 4 + \A py \ 4 ) 



+ y;?1 |A p;c | 2 |A p> ,| 2 + 7p2 (A 2 (A*) 2 + c.c.) 



+K pl (\D x A px \ 2 + \D y A py \ 2 ) + K p2 (\D y A px \ 2 + \D x A py \ 2 ) 
+K p3 ((D x A px T(D y A py ) + c.c.) + K p4 ((D y A px T(D x A py ) + c.c.)}, 
F ST = J J 2 r[r3(|A„,| 2 + lA^.^IA.I 2 + y 4 (\A px \ 2 + |A pv | 2 )|A/ 

+y 5 {(A 2 x + A 2 y )(A:) 2 + c.c.} + 76 {(A 2 px + A 2 y )(A*) 2 + c.c.} 
+ yi {\A px \ 2 - |A p> .| 2 )(A:A d + c.c.) + y 8 {(A 2 A . - Aj^A* + c.c.}}, 
F M = J d 2 r[a m m 2 + fi m m 4 + K m (Vmf], 

F S m = J d 2 r(y ms m 2 \A s \ 2 + y md m 2 \A d \ 2 ), 

F T m = J d 2 ry mp m 2 (\A px \ 2 + \A py \ 2 ), 

F STM = fd 2 r{K spm m[A s ((D x A px T + (D y A py )*) - ((D x A s )A* px + (D y A s )A* py )] 

+K dpm m[A d ((D x A px y - (D y A py y) - ({D x A d )A* px - (D y A d )A* py )] + c.c.}. 
Here, F s , F T , and F M are the free energy for the singlet and triplet SCOPs and the magneti- 
zation, respectively, while F ST , F SM , F TM , and F STM describe their couplings. Note that F is 
invariant under all the symmetry operations of the square lattice. 23) The cubic term F STM has 
derivative couplings of singlet and triplet SCOPs with magnetization, so that triplet (singlet) 
SCOPs would be induced once ferromagnetism coexists with singlet (triplet) superconduc- 
tivity inhomogeneously, as already noted. (In other words, the triplet (singlet) SCOP would 
not be induced if the coexistence occurs uniformly.) This gives a simple and clear inter- 
pretation for previous theoretical results using the BdG method 5) and quasiclassical Green's 
functions, 9 ' 10) in which the occurrence of p-wave SCOPs near the interface between a singlet 
superconductor and a ferromagnet was pointed out. 

Dahl and Sudb0 24) derived the GL free energy from a model with a spin generalized BCS 
term and a Heisenberg exchange term, which is different from ours. They found a cubic term 
in the GL free energy that couples a nonunitary SCOP with magnetization. 
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4. Case of Superconductivity and Antiferromagnetism 

Next, we consider the coexistent states of superconductivity and antiferromagnetism. In 
the AF state the magnetization ntj is oscillating, so we expect that the triplet SCOP will be 
induced even in a uniform AF state once the coexistence occurs. As a slowly varying OP to be 
considered in the continuum theory, we define the staggered magnetization Mj = m ; e iQr/ with 
Q = (n,n). If we assume that the singlet component of SCOP, A (S \ is also slowly varying, 
then the triplet component A <T) should oscillate, as can be seen from eqs. (1 1)-(13). Therefore, 
we define the ^-triplet SCOP A ( * T \j) = A { J\j)e iQ ' r ' . Rewriting eqs. (11)-(13) in terms of M 
and A (7r7) , we find that all terms in these equations do not have staggered oscillations. Defining 
the p x and p y components of the ^--triplet SCOP as 

A£ r) o-)= \{$ T \j)+^m, (21) 

we carry out calculations similar to those in the ferromagnetic case to get GL equations and 
the GL free energy. Here, we present only the resulting expressions for the free energy F AF : 

pAF _ j? , r-AF , pAF , t?AF , r-AF , t?AF , t?AF 

_r S + r T +r ST +r M S M TM S TM' 

Ff = J d 2 r[a pl (KP\ 2 + |A^| 2 ) + & p2 {A%\$?y + c.c) + &(|A,/ + |A,/) 

+ f lA \A% T) \ 2 \A% T) \ 2 + y p2 ((A^) 2 (A<f *) 2 + c.c.) 
+y p3 (\A% T '\ 2 + |A<f | 2 )(A^(A(f )• + c.c.) 
+ K pl {\D x Af\ 2 + \D y A^\ 2 ) + K p2 (\D y AfP\ 2 + \D X A^\ 2 ) 
+K p3 ((D x A$Py(D x A% T) ) + (D y A%PnD y A%P) + c.c.)], 
Ff T = J d 2 r[y 3 (\A% T) \ 2 + |A^ r f)|A,| 2 + y 4 (|Ajf f + |A<f f)|A/ 

+f 5 {((Ajf >) 2 + (A<f ) 2 )(A:) 2 + c.c.} + %{((AfP) 2 + (Ag r) ) 2 )(A*) 2 + c.c.\ 
+HKP\ 2 - |Ag r f)(A:A, + c.c.) + yMtfP? - (Ag r )) 2 )A:A* + c.c.}} 
+y 9 ((tfP)*A% T) + c.c.)\A s \ 2 + y 10 ((Ajf >)*Ag^ + c.c.)\A d \ 2 
+ ril (A^A(f (A:) 2 + c.c.) + y 12 (Ajf >Ajf (A*) 2 + c.c.)], 

F A M F = J d 2 r[a m M 2 + /3 m M 4 + K m (VM) 2 ], 
F sm = J d 2 r(f ms M 2 \A s \ 2 + y md M 2 \A d \ 2 ), 

Ff M = J d 2 v[y mpl M 2 (\A^P\ 2 + \A%p 2 ) + y mp2 M 2 (A% T \A^)* + c.c.)], 
F A s F tm = f d 2 r[y spm MA s (A<«P + + y dpm MA d (AfP - A^f + c.c], 

(22) 
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where F s is the same as that in the ferromagnetic case. The expressions for the coefficients 
appearing in F AF are summarized in Appendix D. The cubic term F A j M in this case couples 
A (S) , A^ T \ and M directly without derivatives. Then the (p x + /? v )-wave [(p x - p y )-vjave] 
^■-triplet component would be induced when s-wave (rf-wave) superconductivity and antifer- 
romagnetism coexist, even in a uniform case. This is consistent with the results of previous 
mean-field calculations that predict the occurrence of the 7r-triplet component in uniformly 
coexistent states of d- wave superconductivity and antiferromagnetism. 2529) This term also 
gives a simple explanation for the occurrence of the ^--triplet SCOP near the interface between 
a singlet superconductor and an antiferromagnet, which was found in numerical calculations 
based on the BdG method. 5) 

5. Summary and Discussion 

We have derived GL equations and the GL free energy for the coexistent states of super- 
conductivity and magnetism microscopically from the extended Hubbard model with on-site 
repulsive and nearest-neighbor attractive interactions. It was found that, in the GL free energy, 
a cubic term that couples singlet and triplet SCOPs with magnetization exists. Owing to this 
term, triplet SCOPs would be induced when ferromagnetism coexists with singlet supercon- 
ductivity inhomogeneously. This gives a simple explanation for previous theoretical studies 
on bilayer systems composed of a ferromagnet and a singlet superconductor. 5 9 ' 10) In the co- 
existent state of antiferromagnetism and singlet superconductivity, ^--triplet SCOPs would be 
induced. This occurs not only in inhomogeneous cases but also in spatially uniform states. 

The validity of the model employed in this paper is limited because of the absence of the 
SU(2) symmetry in spin space. For more general and precise argument of the symmetry of 
the induced OPs, theoretical investigations based on the model that respects this symmetry 
will be necessary, although the present study may capture some of the important aspects. 

In order to study the material dependence of interface states more generally, it is neces- 
sary to derive GL equations from other microscopic models. For example, the low-energy 
electronic states of high- T c cuprate superconductors are described by the t — J model, 30) and 
so the interface state of heterostructures made of high-T c cuprates and magnetic materials 
may be studied using the GL equations derived from this model. 

Numerical study of the GL equations derived from different microscopic models may 
clarify the material dependence of the interface states of heterostructures composed of various 
superconductors and magnetic materials. This problem will be examined separately. 
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Appendix A: Functions Appearing in GL Equations 

The functions L {n) (n = 1, • ■ • ,8) appearing in eqs. (12)-(14) are defined as follows: 

G (k, j + j], -is n )G Q (k + 6, j, ie n ), 

s„ k,S 

L (2) (j, k, k', 7], 6) = VUT 2^ [G (k, j + tj, -is n )G (k + 6, k', is n )G (k', j, is n ) 

-G (k', j + T], -is n )G (k, k', -is n )G (k + 6, j, ie n )\, 
L< 3 \j, k, k', k", 77, 6.6', 6") = -V 3 T ^ G (k, j + J], -is n )G (k + 6, k', is n ) 

xG Q (k", k' +6', -is n )G (k" + 6", j, is n ), 
L^(j, k, k', k", i/, 6) = VU 2 T £ [G (k, j + 77, -is n )G (k + 6, k', is n ) 

xG (k',k", is n )G (k",j, is n ) 

-G (k', j + 77, -is n )G (k, k', -is n )G (k + 6, k", is n )G (k", j, ie n ) 
+G (k", j + 77, -is n )G (k', k", -ie n )G (k, k', -ie n )G (k + 6, j, 
L (5) (j, k) = -UT2_j Go(j, k, ie n )G (k, j, ie n ), 

L®a k, k', 6, 6') = V 2 T G (j, k, is n )G (k', k + 5, -is n )G Q (k' + 6', j, is n ), 

L (7 \j, k, k', k") = -U 3 T 2^ GqU, k, is n )G (k, k', is n )G (k', k", is n )G (k", j, is n ), 

L (8) (j, k, k', k", 6, 6') = V 2 UT ^ [g (j, k, is n )G (k', k + 6, -ie n ) 

xG (k' + 6', k", ie n )G (k", j, ie n ) 

+G (j, k", is n )Go(k", k, is n )G (k', k + 6, -ie n )G (k' + 6', j, is n ) 

-G (j, k, ie n )G (k", k + 6, -is n )G (k', k" , -is n )G (k' + 6', j, ie n j\. 

(A-l) 

Appendix B: Derivation of the Second Order Terms in GL Equations 

In this appendix, we show how to calculate the second-order terms in the GL equations 
for ferromagnetism and superconductivity. Here, the equation for A s is treated as an example. 
(Other OPs can be treated similarly.) The term to be considered is 

7 VUT ^ ^ {A ( p(k))*m k , [<5 (k, j + 77, -is n )G (k + 6, k', is n )G (k', j, is n ) 

v e„ k,k',S (BT) 

-G (k', j + 77, -is n )G (k, k', -is n )G (k + 6, j, «£„)]. 
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We substitute eq. (15) for Af(k) and use a similar approximation to m^, and denote Tj — > r, 
Yk — > r' and i> — > r" . The term without derivatives and A is given by 

l -VUm{v) T^^YuYu ^ E (Af ) (r))^-' p - 5 [Gg(p, * e „)G (p, -fe») 

e„ p x) 6 (B-2) 

-G (p, is„)Go(p, -/£„)], 
where Af is the total number of lattice sites. This term is seen to vanish by putting e n — > -s n 
in the second line. Next the terms that are first order in derivatives and A are given as 

+Af ) (r)(r" - r)^m(r)] — ^ G (pi, -is n )G (p 2 , is n )G (p 3 ,is n ) ( B " 3 ) 

Xg i(Pi+P2)-r' e i(P3-P2)r" g -!(pi+p3)-r/ g -ipr)7 g ip2-5 _ e iPr<5g-iP3^ 

Substituting the relations 

(r' - r)„e i(Pl+P2)r ' = e^^W - ,-_) e «(Pi+B2>fr'-f) > 

<V (B4) 
(r" - r)„e' (P3 - p2) r = e i{P3 - p2) r ( - i—)e^-^< T ~ r \ 

we carry out the integrations over r' and r" after performing the partial integration on p^. 
Then eq. (B-3) becomes 

11 Q 
~1 VU Yj T YjYjnYj [(Af ) (r))*V, 1 m(r)G (p, -w n )G Q (p, w n )i—G Q (p, is n ) 

+m(r)(D ll A'P(T))"G (p, -is,)i-^-Gl(p, fe„)](«T* V"* s - e*'e* s ) 

- 2 E E - ^) 

p "I* <5 

x[(Af>(r))*Vi(r) + 2m(r)(D M Af } (r))*] 

= - 2yc/ ^ E ^(p)^^[ (V ^ m(r)) ( A ^( r )^ + A ;>( r )^y) 

+2ro(r){(D /1 A p ,(r)) V + (D^r)) *w v }]. 

(B-5) 

With the definitions of ^. s/)m and the function h(p) in Appendix C, the last expression is seen 
to give the second-order term appearing in eq. (16). 

Terms that are second order in derivatives and A can be shown to vanish by carrying out 
similar calculations. 
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Appendix C: Coefficients in GL free energy for ferromagnetism and superconductivity 

The coefficients appearing in GL equations [eqs. (16)-(19)] and the GL free energy [eq. 
(20)] are given as follows: 

a m = 4V(\-^Yjh(p)co 2 s(d) ), 

p 

p 

y x = 32 V 4 - I 2 (p)co 2 co 2 d , y 2 = ^yi, 

p 

*p=w{i- 7 ^Yj h< j> )a) % 

p 

fi P = ^V 4 ^J]l 2 (p)co x , 
^p 

7 p i = 32V 4 - ^ h(p)co 2 x co 2 y , y p2 = -y pU 
p 

K P s + K p4 = 4V 2 i £ / 2 (^)(^)(^)^, 
{73, 74, 7i) = 32V 4 — ^ I 2 (p){o) 2 a) 2 x , to d u) 2 x , u) s a) d a) 2 x }, 



p 

1 1 1 

75 = -t73 , 76 = -^74, 7s = -rTz. 

a m = £/(l + ^ 



u 4 



p 



\2N 



p 



*--^Z(|)V"(p>. 



(CI) 



{7m„ 7mrf, 7mp} = ~W 2 U 2 ^- J]^h(P) ~ h(p)){oj], U) 2 , 0J% 

p 

2 1 d^ p 
K spm = 8V U— 2^h(p)-Q^-o) s a) x , 

K dpm = 8V 2 t/-J- V I 3 (p)^L(o d (o x , 
N ^ dp x 
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where a> s = cos p x + cos p y , a> d = cos p x - cos p y , and a> x(y ) = sin p x{y) , and the summa- 
tion on p is taken over the first Brillouin zone. The functions I\, I 2 , and 1$ are defined 
as h(p) = T^G (p,is n )Go(p,-is n ), I 2 (p) = T^^G^(p,ie n )Gl(p,-is n ), and h(p) = 

T ^ G 3 (p, is n )Go(p, -is n ), with f(g p ) being the Fermi distribution function. 

En 
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Appendix D: Coefficients in GL Free Energy for Antiferromagnetism and Supercon- 
ductivity 

The expressions of the coefficients in F AF [eq. (22)] are given as follows: 

a p i = 2V(l --jT^j COs2 Px )' 
p 

4V 2 

a P 2 = — — 2j ^(P) COS Px cos Py, 

p 

P P = SV 4 ^J]l 5 (p)cos 4 p x , 
^p 

%x = 32V 4 - ^ h(p) cos 2 p x cos 2 p y , y p2 = -y p i, 

p 

7 P 3 = 16y4 ^ hip) cos 3 p x cos p y , 
p 

p 

K p3 = -2V — 2_ i I ^P\~) cosp x cosp y , 
p " x 

73 ( 4) = 32V 4 - 2_ J I 6 (p)oj z m cos 2 p x , 

p 

75(6) = 8 V 4 - ^ h(pW s(d) cos 2 p„ 
p 

y 7 = 32V 4 — ^ 7 6 (p)a>,^ cos 2 (D4) 
p 

y 8 = 16V 4 — ^/ 7 (rt^^COS 2 p A -, 

P 

79 (10) = 32 V 4 — h(pW s{d) cos p A - cos p y , 

p 

7ii(i2) = l6y4 Jj Yjh(p)o> 2 s(d) cosp x cosp y , 
p 

s '" = c/ ( 1 + ^Z /8(p) )' 

p 

- U 4 



p 

{fms, Jmd] = -4^ 2 t/ 2 ^ £[2/ 10 (p) + / 7 (p)]{«& O; 2 }, 

P 

{7mpl, 7mp2} = - 4 ^ 2f/2 ^ J] [2^11 (?) + /7(p)]{C0S 2 p. v ,C0Sp. v C0Sp 3 ,}, 

P 

{y V m, 7W = 8V 2 U— Y I n(p) cos p x {co s ,co d }, 

p 
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where the functions appearing in the integrands are defined as 



hip) = 


T ^ G (p, -is n )G Q ip + Q, is n ), 


hip) = 


T 2^ G 2 ip, is n )Gl(p + Q, -is n ), 


hip) = 


T 2^ Glip, ie n )G (p, -is n )G (p + Q, -ie n ), 


hip) = 


T ^ G ip, is n )G Q (p, -is n )G (p + Q, is n )G (p + Q, -is n ), 


hip) = 


T ^ G (p, is n )G ip + Q, i£n), 


hip) = 


T 2^ Glip, ie n )Glip + Q, is n ), 


hoip) = 


T 2j Glip, ie n )G ip, -is n )G Q ip + Q, is n ), 


hlip) = 


T 2^ G 2 ip, is n )G Q (p + Q, ie„)G ip + Q, -is„), 


hlip) = 


T 2^ Goip, is„)G ip, -is n )G Q ip + Q, is n ). 
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